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The goal of hadron physics is the provision of a quantitative explanation of the 
properties of hadrons through a solution of quantum chromodynamics (QCD). This is a 
fundamental problem that is unique in the history of science: never before have we been 
confronted by a theory whose elementary excitations are not those degrees-of-freedom 
readily accessible via experiment. Moreover, QCD generates forces which are so strong 
that less-than 2% of a nucleon's mass can be attributed to the so-called current-quark 
masses that appear in the QCD Lagrangian; viz., forces that generate mass from almost 
nothing, a phenomenon known as dynamical chiral symmetry breaking (DCSB). 

Neither confinement nor DCSB is apparent in QCD's Lagrangian and yet they play 
the dominant role in determining the observable characteristics of real-world QCD. The 
physics of hadrons is ruled by emergent phenomena, such as these, which can only be 
elucidated through the employment of nonperturbative methods in quantum field theory. 
This is both the greatest novelty and the greatest challenge within the Standard Model. 
We must find essentially new ways and means to explain precisely via mathematics the 
observable content of QCD. Herein we provide a snapshot of recent progress in hadron 
physics made using QCD's Dyson-Schwinger equations (DSEs). The complex of DSEs 
is a powerful tool, which has been employed with marked success to study confinement 
and DCSB, and their impact on hadron observables. This is apparent from Ref. [1]. 

Asymptotic coloured states have not been observed, so no solution to QCD will be 
complete if it does not explain confinement. This means confinement in the real world, 
which contains quarks with light current-quark masses. That is distinct from the artificial 
universe of pure-gauge QCD without dynamical quarks, studies of which tend merely to 
focus on achieving an area law for a Wilson loop and hence are irrelevant to the question 
of light-quark confinement. 

Confinement can be related to the analytic properties of QCD's Schwinger functions 
[2] and can therefore be translated into the challenge of charting the infrared behavior 
of QCD's universal /3 -function. This is a well-posed problem whose solution can be 
addressed in any framework enabling the nonperturbative evaluation of renormalisation 
constants. In this connection we note that the hadron spectrum [3], and hadron elastic 
and transition form factors [4, 5] provide unique information about the long-range inter- 
action between light-quarks and, in addition, the distribution of a hadron's characterising 



properties amongst its QCD constituents. However, to make full use of extant and forth- 
coming data, it will be necessary to have Poincare covariant theoretical tools that enable 
the reliable study of hadrons in the mass range l-2GeV. Crucially, on this domain both 
confinement and DCSB are germane; and the DSEs provide such a tool. 

DCSB; namely, the generation of mass from nothing, does take place in QCD. It 
arises primarily because a dense cloud of gluons comes to clothe a low-momentum 
quark. This is best seen by solving the DSE for the dressed-quark propagator [6]; i.e., 
the gap equation. However, the origin of the interaction strength at infrared momenta, 
which guarantees DCSB through the gap equation, is unknown. This relationship ties 
confinement to DCSB. The reality of DCSB means that the Higgs mechanism is largely 
irrelevant to the bulk of normal matter in the universe. Instead the single most important 
mass generating mechanism for light-quark hadrons is the strong interaction effect of 
DCSB; e.g., one can identify it as being responsible for 98% of a proton's mass. 

In chiral-limit QCD, DCSB is most basically expressed in a strongly momentum- 
dependent dressed-quark mass; viz., M{p^) in the quark propagator: 

Si,) = I = ^(P') (1) 

iY-pA{p^)+B{p^) iY-p+M{p^)' 

The appearance and behaviour of M{p^) are essentially quantum field theoretic effects, 
unrealisable in quantum mechanics. The running mass connects the infrared and ultravi- 
olet regimes of the theory, and establishes that the constituent-quark and current-quark 
masses are simply two connected points on a single curve separated by a large momen- 
tum interval. QCD's dressed-quark behaves as a constituent-quark, a current-quark, or 
something in between, depending on the momentum of the probe which explores the 
bound-state containing the dressed-quark. It follows that calculations addressing mo- 
mentum transfers > M^, where M is the mass of the hadron involved, require a 
Poincare-covariant approach that can veraciously realise quantum field theoretical ef- 
fects [4]. Owing to the vector-exchange character of QCD, covariance also guarantees 
the existence of nonzero quark orbital angular momentum in a hadron's rest-frame [7,8]. 

Through the gap and Bethe-Salpeter equations (BSEs) the pointwise behaviour of 
QCD's /3-function determines the pattern of chiral symmetry breaking. Since these and 
other DSEs connect the /3-function to experimental observables, the comparison be- 
tween computations and observations of the hadron properties can be used to chart 
the /3 -function's long-range behaviour. In order to realise this goal a nonperturbative 
symmetry-preserving DSE truncation is necessary. Steady progress has been made with 
a scheme that is systematically improvable [9, 10]. On the other hand, significant quali- 
tative advances in understanding the essence of QCD could be made with symmetry- 
preserving kernel Ansatze that express important additional nonperturbative effects, 
which are impossible to capture in any finite sum of contributions. 

Such an approach is now available [11]. It begins with a novel form for the axial- 
vector BSE, which is valid when the quark-gluon vertex is fully dressed. Therefrom, a 
Ward-Takahashi identity for the Bethe-Salpeter kernel is derived and solved for a class of 
dressed quark-gluon-vertex models. The solution yields a symmetry -preserving closed 
system of gap and vertex equations. As the analysis can readily be extended to the vector 
equation, a comparison is possible between the responses of pseudoscalar- and scalar 



TABLE 1. Axial-vector and vector meson masses calculated in three truncations of the coupled 
gap and Bethe-Salpeter equations. The last column was obtained using the standard Ball-Chiu An^afz 
augmented by the quark anomalous chromomagnetic moment in Eqs. (3), (4). 
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meson masses to nonperturbatively dressing the quark-gluon vertex. The resuh indicates 
that spin-orbit splitting in the meson spectrum is enormously enhanced by DCSB. 

It has been conjectured [11] that the full realisation of DCSB in the Bethe-Salpeter 
kernel will have a material impact on mesons with mass greater than 1 GeV. Moreover, 
that it can overcome a longstanding failure of theoretical hadron physics. Namely, no 
extant hadron spectrum calculation is believable because all symmetry preserving stud- 
ies produce a splitting between vector and axial-vector mesons that is far too small: 
just one-quarter of the experimental value (see, e.g., Refs. [12, 13, 14]). In this connec- 
tion, preliminary results are now available [1] and they are listed in Table 1. The second 
numerical column reports results obtained in rainbow-ladder truncation; viz., leading- 
order in the systematic and symmetry-preserving truncation scheme of Ref. [10]. As 
anticipated, while the p-meson mass is acceptable, the ai-mass is far too small. 

The procedure introduced in Ref. [11] enables meson masses to be calculated using 
a symmetry-preserving DSE truncation whose diagrammatic content is unknown. One 
can therefore elucidate the effect of an essentially nonperturbative DCSB component 
in dressed-quark gluon vertex on the p-ai complex; in this case, the Ag term in the 
Ball-Chiu vertex [15], which had an enormous impact in the scalar channel: 

ir^{q,k) = ii:A{q\k^)rn+2in[ir-i^A{q\k^)+AB{q\k^)], (2) 

where i = {q + k)/2, L^{q,k-P) = [^{q-P)+^{k;P)]/2 and A^{q,k;P) = Mq;P) - 
^{k; P)]/ [q^ — k^] . The results obtained with this vertex are shown in Table 1 : the DCSB 
Ag-term boosts the <3i mass, which is a positive outcome, but it simultaneously boosts the 
p mass, such that the mass- splitting is practically unchanged from the rainbow-ladder 
result. Was Ref. [11] too optimistic in imagining that the new scheme could provide the 
first realistic meson spectrum encompassing states with mass greater-than 1 GeV? 

Before answering, let us return to a consideration of chirally symmetric QCD. That 
theory exhibits helicity conservation so that, perturbatively, the quark-gluon vertex can- 
not have a term with the helicity-flipping characteristics of Ag. There is another feature 
of massless fermions in gauge field theories; namely, they cannot posses an anoma- 
lous chromo/electro-magnetic moment because the term that describes it couples left- 
and right-handed fermions. However, if the theory's chiral symmetry is strongly bro- 
ken dynamically, why shouldn't the fermions posses a large anomalous chromo/electro- 
magnetic moment? Such an effect is expressed in the quark-gluon-vertex via a term 

T^''^{q,k) = Oi,y{q-k)yX5{q,k) (3) 



where, owing to DCSB, the natural strength is represented by the Ansatz 

T5{q^)=ABiq\k^). (4) 

NB. Based on the models in Refs. [4, 16], 2M£Z(M|)Afi(M|,M|) ~ -i, where Me is 
the Euclidean constituent-quark mass, defined in Ref. [17]. 

Using the procedure introduced in Ref. [11], the vector and axial- vector vertex equa- 
tions can be solved using the dressed-quark-gluon vertex obtained as the sum of Eqs. (2) 
and (3). The effect is remarkable: the anomalous chromomagnetic moment leads to ad- 
ditional repulsion in the ai channel but significant attraction in the p channel such that, 
for the first time, a realistic result is simultaneously obtained for the masses in both these 
channels, and hence the ai-p mass-splitting - see Table 1 . Furthermore, the origin of the 
splitting, in an interference between As in Eq. (2) and T5 in Eqs. (3), (4) is intuitively 
appealing. In the chiral limit the mass- splitting between parity partners should owe 
solely to DCSB and here that is seen explicitly: in the absence of DCSB, Ag = = T5. 
The rainbow-ladder result is also understood: this truncation fails to adequately express 
DCSB in the Bethe-Salpeter kernel and hence cannot realistically split parity partners. 

Table 1 is a "first-guess" result; i.e., there was no tuning of the strength in Eq. (4), 
so how reliable can it be? This question amounts to deciding whether a realistic size 
is assumed for a light-quark's anomalous chromomagnetic moment. Fortunately, an 
analysis is available of results for the dressed-quark-gluon vertex obtained through 
numerical simulations of quenched-QCD [18]. This study shows that T5 is dynamically 
two orders-of-magnitude larger than the one-loop perturbative result and, indeed, is of 
the same magnitude and possesses the same domain of significant support as Ae(^^, k^), 
precisely in accordance with the assumption we have made. In addition, efficacious 
model studies support a moment of similar size [19, 20]. 

At this point it is natural to consider whether DCSB in QCD can also generate a 
large quark anomalous electromagnetic moment term in the quark-photon vertex. In 
perturbation theory, of course, since it doesn't express DCSB, the quark's anomalous 
electromagnetic moment is small [21]. One obtains the same answer in the rainbow- 
ladder truncation; e.g., the F(, and terms in Ref. [22], which combine to form T5 in our 
notation, contribute less-than 1% to the pion's electromagnetic form factor. However, 
as we've already seen, this truncation doesn't adequately incorporate DCSB into the 
Bethe-Salpeter kernel. Only the method of Ref. [11] can readily provide an answer to 
this question and we are currently compiling the necessary information. 

Another important problem is the computation of the parton distribution functions 
of the most accessible hadrons [23]. In connection with uncovering the essence of the 
strong interaction, the behaviour of the valence-quark distribution functions at large 
Bjorken-x is most relevant. Owing to the dichotomous nature of Goldstone bosons, 
understanding the valence-quark distribution functions in the pion and kaon is of great 
importance. Moreover, given the large value of the ratio of s-to-u current-quark masses, 
a comparison between the pion and kaon structure functions offers the chance to chart 
effects of explicit chiral symmetry breaking on the structure of would-be Goldstone 
modes. There is also the prediction [24, 25] that a theory in which the quarks interact 
via \/k^ vector-boson exchange will produce valence-quark distribution functions for 
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FIGURE 1. Left panel. Valence M-quark distribution functions, computed as described in the text: solid 
curve, Mv(x); and dashed curve, u^{x). Applying leading-order QCD evolution from Qq — 0.32GeV^ 
to = 25GeV^, explained in Sec.II.D of Ref. [23], one obtains the other two curves from these 
starting distributions: dashed curve, u^{x;25); and dotted curve, uf{x,25). Right panel. Model ratio 
Uy /uy evaluated at Qq = 0.32GeV^, Dashed curve, and = 25GeV^, solid curve, /m" = 0.13. 
Under the right conditions, u^ /u^ should equal the ratio of kaon-to-pion Drell-Yan cross-sections, and 
we reproduce that obtained from a sample of dimuon events with invariant mass 4.1 < M < 8 .5 GeV [26] . 



which 

^vWoc(l-x)2+r,x^0.85, (5) 

where 7 ^ is an anomalous dimension that grows with increasing momentum transfer. 
(See Sec.VLB.3 of Ref. [23] for a detailed discussion.) 

Experimental knowledge of the parton structure of the pion and kaon arises primarily 
from pionic or kaonic Drell-Yan scattering from nucleons in heavy nuclei [26, 27]. 
Theoretically, given that DCSB plays a crucial role in connection with pseudoscalar 
mesons, one must employ an approach that realistically expresses this phenomenon. 
The DSEs therefore provide a natural framework: a study of the pion exists [28] and one 
of the kaon is underway [29]. 

The results to be anticipated from the latter study have been illustrated [1] through an 
internally consistent calculation based upon the QCD-improvement of a simple model 
used already for pion and kaon distribution functions [30]. In Fig. 1 we depict our 
computed distributions themselves and relevant ratios. Aspects of the curves are model- 
independent. For example, owing to its larger mass, one anticipates that the 5-quark 
should carry more of the charged-kaon's momentum than the M-quark. This explains 
why the support of xu^{x) is shifted to lower-x than that xu^{x). QCD evolution is an 
area-conserving operation on the distribution function, which shifts support from large-jc 
to small-x. Thus, while both Uv''^{x;Qo) °^ (1 —x)^ for x ~ 1, 

K-''{x;Q)"^\l-xr,a = 2J. (6) 

These observations explain the qualitative behaviour of the evolved distributions. Con- 
cerning the ratio, as a consequence of the form of the evolution equations its value at 



X = 1 is invariant under evolution, whereas the value at x = approaches one under 
evolution owing to the increasingly large population of sea-quarks produced thereby. 

There are many reasons why this is an exciting time in hadron physics. We have 
focused on one. Namely, through the DSEs, we are positioned to unify phenomena as 
apparently diverse as the: hadron spectrum; hadron elastic and transition form factors, 
from small- to large-Q^; and parton distribution functions. The key is an understanding 
of both the fundamental origin of nuclear mass and the far-reaching consequences of the 
mechanism responsible; namely, DCSB. These things might lead us to an explanation 
of confinement, the phenomenon that makes nonperturbative QCD the most interesting 
piece of the Standard Model. 
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